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Abstract 


This  thesis  investigates  the  use  of  an  improved  finite  element 
model  of  a  T-38  horizontal  stabilizer  for  flutter  analysis  using 
NASTRAN.  The  procedure  for  evaluating  the  effect  of  repairs  on  the 
flutter  speed  is  developed  and  its  sensitivity  to  several  modeling  as¬ 
sumptions  and  practices  is  presented.  The  procedure  is  to  be  used  by 
Air  Force  engineers  to  evaluate  repair  limits  of  T-38  stabs. 

The  NASTRAN  flutter  speed  calculations  are  based  upon  a  structural 
model,  an  unsteady  aerodynamic  model  and  an  interface  model  (splines). 

The  structural  model  consists  of  a  flat  array  of  plates  and  bars  rep¬ 
resenting  the  single  spar  honeycomb  structure.  A  flat  array  of  doublet 
lattice  panels  were  chosen  for  the  unsteady  aerodynamic  model.  An 
interface  between  these  two  models  is  required  to  find  the  dO'«'n  washes 
on  the  aerodynamic  model  due  to  motion  of  the  structural  model.  This 
is  accomplished  by  interpolation  via  surface  splines. 

Solution  sensitivity  to  various  parameters  were  considered  to 
determine  their  effect  on  the  flutter  conditions.  In  addition  to  using 
the  doublet  lattice  method,  strip  theory  and  mach  box  methods  were  used 
in  the  flutter  analysis.  This  allowed  a  comparison  of  flutter  conditions 
versus  aerodynamic  theory  to  be  made. 
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INVESTIGATION  OF  AN  IMPROVED  FINITE  ELEMENT 
MODEL  FOR  A  REPAIRED  T-38  HORIZONTAL 
STABILIZER  FLUTTER  ANALYSIS  USING  NASTRAN 


I.  INTRODUCTION 


Background 


San  Antonio  Air  Logistics  Center  (SAALC)  has  the  primary  respon¬ 
sibility  for  all  engineering  and  maintenance  for  the  Northrop/United 
States  Air  Force  T-38  Talon  supersonic  jet  trainer.  This  responsibility 
includes  determining  whether  or  not  a  damaged  T-38  stabilizer  (stab) 
can  be  repaired  within  the  limits  stated  in  the  T.O.  IT-38A-3.  Be¬ 
cause  the  stab  (Fig.  1)  is  a  flutter  critical  component  on  the  aircraft, 
the  repair  limits  are  established  based  >n  changes  to  the  flutter  speed. 
SAALC  currently  uses  a  method  that  conservatively  predicts  flutter 
speeds.  To  improve  their  flutter  speed  predictions,  SAALC  began  an 
investigation,  in  1975,  for  a  more  accurate  analysis.  The  investigation 
was  subsequently  passed  on  to  the  Air  Force  Institute  of  Technology 
(AFIT)  master  of  science  students. 

In  1979,  Lassiter  (Ref  5)  began  his  investigation  into  the  initial 
development  of  a  structural  model  for  use  in  flutter  analysis  using 
NASA  Structural  Analysis  (NASTRAN).  Lassiter's  structural  model  con¬ 
sisted  of  two  dimensional  quadrilateral  plate  and  bar  elements.  The 
finite  element  structural  model  was  found  to  be  lacking  in  stiffness. 


The  weight,  center  of  gravity,  and  mass  distribution  of  the  model 
agreed  fairly  well  with  that  of  the  actual  stab.  Lassiter  conducted 
rigid  and  flexible  root  analysis  of  the  model  and  compared  the  results 
with  calculated  data  from  Northrop  Aircraft  Incorporation  (NAI).  The 
results  from  the  modal  analysis  showed  good  representations  of  the 
mode  shapes  but  some  of  the  frequencies  varied  as  much  as  20  cycles  per 
second  (cps).  Lassiter  used  the  finite  element  model  with  the  above- 
mentioned  discrepancies  to  perform  a  flutter  analysis  which  yielded  a 
flutter  speed  that  was  well  below  the  NAI-reported  flutter  speed.  The 
NAI  57-59  report  stated  that  the  series  two  and  series  three  stabs  have 
the  same  stiffness  and  Lassiter's  model  of  a  series  three  stab  was  found 
lacking  stiffness.  Since  most  of  the  previous  analyses  and  tests  by  NAI 
were  done  on  a  series  two  stab,  further  investigation  was  necessary  on 
the  series  three  stab. 

Thomson  (Ref  18)  worked  to  verify  the  series  three  finite  element 
model  under  static  loads,  performed  vibration  tests  on  a  stab,  accom¬ 
plished  a  NASTRAN  modal  analysis  on  a  series  three  model,  and  performed 
a  NASTRAN  flutter  analysis.  The  series  two  and  three  stabs  were  modeled 
for  static  load  comparison  and  both  lacked  stiffness  when  compared  with 
NAI's  experimental  results.  The  series  three  model  lacked  more  in 
stiffness,  causing  Thomson  to  increase  the  elastic  and  shear  moduli  by 
30  percent.  The  calculated  frequencies  between  the  series  two  and  three 
models  were  very  close  to  each  other  and  the  mode  shapes  were  good  rep¬ 
resentations  of  his  experimental  mode  shapes. 


For  verification  of  a  steady  aerodynamic  model,  Thomson  made  a 
comparison  between  USSAERO  and  NASTRAN  models,  and  preliminary  exper¬ 
imental  data.  There  was  good  correlation  between  the  models  and  the 
experimental  data.  Thomson  concluded  that  NASTRAN  modeled  the  aerody¬ 
namics  of  the  airfoil  correctly  for  the  steady  case.  The  results  of 
his  flutter  analysis  was  well  below  reported  flutter  speeds.  Thomson 
recommended  that  corrections  to  the  model  (tuning)  should  be  accom¬ 
plished  and  the  torsional  spring  stiffness  of  the  hydraulic  actuator 
system  should  be  investigated. 

Sawdy  (Ref  17)  used  Thomson's  two-dimensional  model  for  the  pre¬ 
diction  of  static  displacement  for  the  T-38  stab.  He  completed  exper¬ 
imental  static  deflection  tests  on  the  series  three  stab  and  performed 
a  static  analysis  on  the  finite  element  model.  The  experimental  test 
set-up  for  the  stab  measured  the  displacements  for  ten  load  conditions 
with  the  displacements  measured  at  twenty-five  locations  for  each 
loading  condition.  The  test  paralleled  that  developed  by  Northrop  for 
the  series  two  stab  with  deviations  in  pad  size,  several  load  conditions, 
the  use  of  a  load  cell  in  place  of  the  control  system  actuator  ram,  and 
boundary  conditions.  Because  of  these  differences,  a  good  comparison  of 
the  static  analysis  of  the  series  two  and  series  three  stabs  could  not 
be  accomplished.  He  explained  the  use  of  optimization  theory  for  tuning 
the  finite  element  model  based  upon  the  static  displacements  found  dur¬ 
ing  his  experimental  tests.  A  value  of  1.18  x  10®  in-lb/rad  for  model¬ 
ing  the  load  cell  was  found  by  direct  search  to  produce  the  lowest 


error  value.  With  the  combined  optimum  torque  tube  bending  Inertia  and 
actuator  arm  spring  constant,  the  error  was  found  to  be  1.6  compared  to 
10.4  before  any  optimization  attempts.  Sawdy  concluded  that  the  finite 
element  model  using  orthotropic  plate  elements  represents  a  good  stab¬ 
ilizer  finite  element  model. 

Concurrent  with  Sawdy 's  investigation.  Dodge  (Ref  2)  began  improv¬ 
ing  the  structural  model's  frequency  response.  Flutter  analysis  re¬ 
quires  accurate  mode  shapes  and  frequencies  from  the  model.  This  is 
accomplished  by  adjusting  (tuning)  the  model  to  match  the  characteris¬ 
tics  of  the  actual  stab.  The  modulus  of  elasticity  and  shear  modulus 
were  the  material  properties  chosen  by  Dodge  for  tuning  the  model .  The 
cross-sectional  properties  chosen  were  the  bending  and  polar  moments  of 
inertia  of  the  spar,  trailing  edge  and  leading  edge  elements.  Changing 
the  structural  mass  and  varying  the  airfoil  thickness  were  also  chosen 
for  a  total  of  six  parameters  used  in  tuning  the  model.  Dodge  concluded 
that  an  increase  of  37  percent  in  airfoil  thickness  was  considered  the 
best  tuned  model  without  changing  the  type  or  number  of  elements  in  the 
model.  The  airfoil  thickness  was  found  to  affect  the  torsional  and 
bending  frequencies  equally.  Modeling  of  the  control  system  pitch 
stiffness  was  found  to  affect  the  first  torsion  mode  the  most.  This 
mode  is  considered  by  Northrop  to  be  the  most  crucial  mode  because  of 
its  effect  on  the  flutter  speed.  Dodge  found  a  pitch  spring  constant 


of  4.40  X  10®  in-lbs/rad  was  necessary  to  tune  the  model  to  NAI's  first 
torsion  mode.  His  best  model  was  the  result  of  increasing  the  airfoil 
thickness  by  37  percent  and  using  a  control  system  spring  constant  of 
4.40  X  10®  in-lbs/rad. 

Statement  of  the  Problem 

Develop  a  procedure  to  select  the  aerodynamic  model,  structural 
model,  and  splines  for  use  in  the  flutter  analysis  process  with  the 
goal  of  determining  flutter  speed  changes  of  a  repaired  stab.  This 
procedure  is  to  be  used  by  Air  Force  engineers  for  updating  current 
repair  limitations  to  the  stab. 


Solution  Approach 


A  normal  modes  analysis  was  accomplished  on  Sawdy's  and  Dodge's 
best  model  in  order  to  determine  the  most  appropriate  model  to  use  in 
the  investigation.  The  structural  model  chosen  was  the  model  used  by 
Dodge  (Fig.  2).  A  complete  modal  analysis  of  the  model  was  accomplished 
to  find  the  frequencies  of  vibration  and  mode  shapes  in  the  frequency 
region  of  interest.  The  aerodynamic  model  was  chosen  based  on  the  re¬ 
quirements  of  the  doublet  lattice  method.  According  to  Huttsell  (Ref 
4),  Noll  (Ref  14),  and  Pollack  (Ref  15),  the  doublet  lattice  method  is 
considered  the  least  conservative  of  all  subsonic  aerodynamic  theories 
and  it  is  used  as  a  standard  in  industry. 


Eastep  (Ref  3)  states  that  the  doublet  lattice  method  is  preferred 
over  other  available  methods  and  it  is  also  considered  less  conserva¬ 
tive  when  compared  with  the  supersonic  theories  in  the  low  supersonic 
range.  If  the  aerodynamic  theory  is  changed,  the  aerodynamic  model 
must  be  changed.  The  required  interface  between  the  structural  and 
aerodynamic  models  is  the  spline.  The  spline  is  the  model  chosen  to 
find  the  down  washes  at  the  intersection  of  the  aerodynamic  panels 
(aero  points)  due  to  motion  at  the  structural  grid  points  and  the 
spline  chosen  is  dependent  upon  the  aerodynamic  theory  used. 

Using  the  above  described  models,  the  NASTRAN  flutter  procedure 
was  used  to  determine  the  flutter  speed  and  frequency.  The  flutter 
conditions  obtained  only  approximate  to  the  actual  flutter  conditions. 
The  actual  flutter  conditions  are  obtained  using  an  iterative  procedure 
called  match  pointing.  The  actual  flutter  speed  and  frequency  are 
obtained  by  varying  the  assumed  mach  number  and  reduced  frequency  until 
the  mach  number  at  flutter  matches  the  assumed  mach  number.  Huttsell, 
Noll,  and  Pollack  state  that  the  general  practice  is  to  match  point  if 
the  predicted  flutter  conditions  are  close  to  the  operational  limits  of 
the  aircraft,  including  a  15  percent  safety  factor.  If  the  flutter 
conditions  are  not  close  to  these  limits,  match  pointing  is  not 
accomplished.  This  practice  is  to  be  followed  during  this  investiga- 


After  determining  the  flutter  conditions,  various  input  parameters 
were  changed  to  determine  changes  in  the  predicted  flutter  conditions. 
Once  this  information  was  obtained,  mass  was  added  to  the  structural 
model  to  simulate  a  repaired  stabilizer  and  the  flutter  procedure  was 
repeated.  In  addition  to  these  procedures,  strip  theory  and  the  mach 
box  methods  were  used  to  obtain  a  comparison  between  aerodynamic 
theories  and  flutter  conditions.  The  aerodynamic  model  was  changed 
and  the  process  repeated  for  these  theories.  Simulated  repairs  were 
accomplished  by  changing  the  structural  model  and  the  flutter  pro¬ 
cedure  repeated  to  determine  changed  flutter  conditions  for  the  stab 
using  both  theories. 


II.  NORMAL  MODES  ANALYSIS 


Modal  analysis  is  a  necessary  step  that  must  be  performed  before 
attempting  a  flutter  analysis.  The  frequencies  and  mode  shapes  of 
the  NASTRAN  structural  model  were  calculated  and  compared  with  the 
results  obtained  by  NAI  calculations  (Ref  9)  and  NAI's  experimental 
data  (Ref  8)  of  the  stab  using  the  flexible  root  installed  boundary 
conditions.  Also,  modal  analysis  is  performed  to  determine  the  modes 
eliminated  and  modified  by  applying  Guyan  reduction.  The  remaining 
frequencies  and  mode  shapes  are  then  used  by  NASTRAN  in  the  flutter 
analysis  procedure.  A  brief  discussion  of  NASTRAN 's  modal  analysis 
is  necessary. 

Rigid  Format  Three  (Ref  13)  and  the  displacement  approach  is  used 
for  the  normal  modes  analysis.  The  four  methods  of  algebraic  eigen¬ 
value  extraction  available  using  NASTRAN  are  the  tridiagonal,  deter¬ 
minant,  inverse  power,  and  the  Upper  Hessenberg.  The  inverse  power 
method  was  chosen  because  the  eigenvalues  and  eigenvectors  could  be 
found  at  the  same  time.  As  a  tracking  method,  the  inverse  power  method 
extracts  the  roots,  one  at  a  time,  by  iterative  procedures  applied  to 
the  original  dynamic  matrices.  For  the  inverse  power  method,  the 
eigenvalue  problem  is  stated  as 


[K  -  xM]  {U}  =  0 


where  K  Is  a  stiffness  matrix,  M  Is  a  mass  matrix,  and  U  Is  a  displace¬ 
ment  vector,  and  X  Is  the  square  of  the  natural  frequency.  With  the 
aid  of  GCSNAST  (Ref  6),  the  mode  shapes  can  be  depicted  graphically. 

For  all  present  analyses,  the  NASTRAN  structural  model  used  Is 
the  best  model  by  Dodge,  unless  otherwise  stated.  A  normal  modes 
analysis  was  performed  and  the  results  were  compared  to  NAI  calculated 
and  experimental  data  and  results  of  Eglin's  Ground  Vibration  Tests 
(GVT)  (Ref  16).  These  results,  based  on  the  Installed  boundary 
conditions,  are  listed  In  Table  I. 


Table  I 

Comparison  of  the  First  Three  Vibration  Modes 
Eglin's  GVT  NAI  Computed  NAI  Experimental  NASTRAN 


Mode 

Freq  (cps) 

Freq  (cps) 

Freq  (cps) 

Freq  (cps) 

1st  Bending 

18.52 

17.61 

17.3 

17.37 

1st  Torsion 

50.20 

44.89 

44.9 

44.33 

2nd  Bending 

70.69 

78.76 

71.7 

77.64 

The  results  from  NASTRAN  compare  favorably  with  the  NAI  calculated  and 
experimental  data.  These  results  Imply  that  the  structural  model  has 


been  tuned  to  match  the  characteristics  of  the  actual  stab. 


Knowing  what  modes  are  used  In  the  flutter  analysis  process  Is 
very  Important  to  ensure  no  unintended  modes  are  used.  This  check  Is 
best  accomplished  using  the  normal  modes  analysis.  The  frequency  range 
of  Interest  Is  given  as  Input  data  and  the  resulting  frequencies  and 
their  order  of  extraction  are  listed.  The  Importance  of  the  extraction 
order  Is  best  presented  by  an  example.  If  three  of  the  frequencies  were 
119,  124,  and  170  cycles  per  second  and  their  extraction  order  was  124, 
119,  and  170,  the  119  cps  would  not  be  Included  In  the  flutter  analysis 
If  the  upper  limit  was  100.  This  Is  because  NASTRAN  will  find  all 
values  In  the  range,  0  to  100,  and  one  value  out  of  that  range  (124) 
which  Is  calculated  first.  The  frequency  of  119  cps  could  be  Included 
In  the  flutter  analysis  and  would  appear  as  an  unknown  mode  to  the 
user. 

NASTRAN 's  flutter  analysis  procedure  does  not  allow  the  use  of 
all  degrees  of  freedom.  The  unwanted  degrees  of  freedom  are  eliminated 
by  using  the  Guyan  reduction  procedure.  The  lateral  displacement  was 
the  only  degree  of  freedom  retained  for  flutter  computations.  The 
normal  modes  analysis  was  repeated  with  and  without  application  of 
Guyan  reduction  to  determine  the  extraction  order,  modes  eliminated, 
and  the  frequencies  that  were  changed.  These  results  are  listed  in 
Tables  II  and  III  and  shown  In  Figures  3  through  9. 
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The  drag  mode  (Fig.  4)  was  eliminated  using  Guyan  reduction,  but 
the  other  frequencies  show  no  significant  change.  This  mode  is  re¬ 
ferred  to  as  the  drag  mode  because  it  is  the  deformation  that  is 
parallel  to  the  wind  direction. 

The  first  torsion  mode  was  reported  by  NAI  (Ref  9)  to  be  the 
flutter  critical  mode  and  it  is  the  mode  that  is  directly  affected  by 
changes  in  the  control  system  pitch  stiffness.  Variations  in  control 
system  pitch  stiffness  were  accomplished  to  determine  the  first 
torsional  modes  sensitivity  to  this  parameter.  Two  of  the  values  used 
were  1.72  x  10®  in-lbs/rad  based  upon  bench  tests  of  the  control  system 
actuator  by  NAI  (Ref  10)  and  4.4  x  10®  in-lbs/rad  from  the  model 
tuning  by  Dodge.  Unless  otherwise  specified.  Dodge's  value  of  4.4  x 
10®  in-lbs/rad  is  used  for  all  flutter  analyses.  Table  IV  lists  the 
changes  in  the  first  torsion  frequency  due  to  changes  in  the  modeled 
control  system  pitch  stiffness. 


Table  IV 

1st  Torsion  Frequency  Change  Due  to  Spring  Constant  Change 


Spring  Constant  (10®) 
1.00 
1.18 
1.72 
2,00 
3.00 


Frequency  (cps) 
29.36 
31.08 
35.11 
36.73 
40.87 


4.40 


44.33 


III.  FLUTTER  ANALYSIS  USING  DOUBLET  LATTICE  METHOD 


In  addition  to  the  structural  model,  one  of  the  most  important 
factors  in  flutter  analysis  is  the  aerodynamic  model.  The  attributes 
of  a  good  aerodynamic  model  depend  on  the  aerodynamic  theory  used. 

With  NASTRAN,  there  are  five  aerodynamic  theories  available;  all  as¬ 
sume  small  amplitude  sinusoidal  motions  and  the  transient  aerodynamic 
forces  are  obtained  by  fourier  methods.  The  aerodynamic  analysis  is 
based  on  the  finite  element  approach  in  the  same  manner  as  structural 
analysis.  The  finite  elements  are  strips  or  boxes  for  which  there  are 
aerodynamic  forces.  The  lattice  methods  utilize  arrays  of  trapezoidal 
boxes  whose  sides  are  parallel  to  the  airflow. 

According  to  Ref  12,  the  aspect  ratio  of  the  trapezoidal  boxes 
should  be  less  than  0.08  times  the  velocity  divided  by  the  greatest 
frequency  of  interest,  but  no  less  than  four  boxes  per  chord  should  be 
used.  This  led  to  an  aerodynamic  model  consisting  of  eight  chordwise 
and  eight  spanwise  boxes  (Fig.  10).  The  stab  is,  therefore,  modeled 
as  a  flat  doublet  panel,  divided  into  an  eight  by  eight  array  of  boxes. 
With  the  structural  and  aerodynamic  models  chosen,  the  interface  model 
can  be  selected. 

Choosing  the  proper  interface  model  is  another  important  factor 
in  the  flutter  analysis  process.  This  interface  is  based  upon  the 
theory  of  splines  (Ref  12).  High  aspect  ratio  or  beam  models  require 


the  use  of  linear  splines.  Since  the  stabilizer  is  modeled  as  a  flat 
plate,  the  use  of  surface  splines  was  appropriate.  Splines  are  used 
to  find  the  down  washes  at  the  aero  points  due  to  motion  at  the 
structural  points.  All  aero  elements  must  be  referred  to  by  a  spline 
or  that  element  will  have  no  wash.  If  desired,  the  spline  can  be 
segmented  into  any  rectangular  subarray  of  boxes  on  a  panel.  The 
input  points  to  the  spline  are  the  desired  structural  grid  points  to 
be  used  in  the  analysis.  With  the  spline  defined,  all  models  required 
for  flutter  analysis  have  been  specified.  These  models,  structural, 
aerodynamic,  and  splines,  will  not  be  changed  in  the  analyses  unless 
specifically  stated. 

Before  the  flutter  analysis  process  can  be  used,  the  method  of 
modal  flutter  analysis  must  be  chosen.  Modal  flutter  analysis  methods 
available  using  NASTRAN  are  the  K,  KE,  and  PK.  The  K  method  allows 
looping  through  three  sets  of  parameters:  density  ratio,  mach  number, 
and  reduced  frequency.  The  KE  method,  though  more  efficient  than  the 
K  method,  is  restricted  to  no  damping  matrix  and  no  eigenvector  recov¬ 
ery.  The  PK  method  tracks  the  aerodynamic  matrices  as  frequency  de¬ 
pendent  springs  and  dampers,  where  the  frequency  is  estimated  by  the 
user  and  the  eigenvalues  are  found.  The  K  method  was  chosen  because 
eigenvector  recovery  is  available  for  mode  shape  description  at  the 
flutter  condition.  To  determine  the  flutter  condition,  a  plot  of 
velocity  (V)  versus  damping  (g)  for  various  reduced  frequencies  (k) 
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is  accomplished  to  locate  the  point  where  the  damping  is  zero.  The 
velocity  at  this  point  is  the  flutter  velocity  (V^)  and  the  frequency 
is  the  flutter  frequency 

Using  the  K  method,  the  tuned  structural  model,  the  eight  by  eight 
aerodynamic  model,  and  a  surface  spline  with  inputs  from  all  structural 
nodes  (Fig.  11),  the  flutter  analysis  process  was  accomplished.  All 
analyses  are  based  upon  sea  level  conditions  and  .5  Mach  unless 
stated  otherwise.  The  flutter  condition  from  the  analysis  was  a  of 
916  Knots  and  of  29.7  cps.  These  results  compared  favorably  with 
the  NAI  (Ref  11)  reported  of  830  Knots  and  of  29  cps.  The  report 
also  stated  that  the  830  knot  flutter  speed  was  ten  to  twenty  percent 
below  expected  values.  A  V-g  diagram,  at  various  values  of  reduced 
frequency,  was  plotted  (Fig.  12)  for  the  first  three  modes  and  compared 
to  a  similar  plot  by  NAI  (Fig.  13).  Additionally,  modes  4  through  6 
were  plotted  on  a  V-g  diagram  (Fig.  14)  as  well  as  a  plot  for  the 
first  six  frequencies  versus  velocity  (Fig.  15).  Figure  12  was  almost 
identical  to  those  produced  by  NAI  but  with  the  noted  difference  in 
flutter  speed.  NAI's  830  Knot  was  based  on  three  percent  damping 
whereas  the  results  in  this  report  are  based  on  zero  damping. 

Changes  to  these  models  and  other  input  variables  were  investigat¬ 
ed  to  determine  their  effect  on  V^.  These  changes  were  varying  the 
modeled  control  system  pitch  stiffness,  segmenting  the  splines. 


re-'ucing  the  number  of  aerodynamic  boxes,  and  varying  the  altitudes. 

The  previously  mentioned  values  for  the  control  system  pitch  stiffness 
were  used  in  the  flutter  analysis  procedure  to  determine  the  effect  on 
V^.  The  results,  listed  in  Table  V,  confirm  the  expected  flutter  speed 
changes  due  to  a  change  in  the  first  torsion  frequency.  The  change  in 
the  first  torsion  frequency  is  directly  effected  by  the  control  system 
pitch  stiffness  resulting  in  a  flutter  speed  decrease  as  the  torsional 
frequency  decreases. 

Table  V 

Changes  Due  to  a  Spring  Constant  Change 

In  Vacuum 

Spring  First  Torsion  „ 

Constant  Freq  (cps)  *f  (knots)  ‘^f  (cps) 


1.00 

29.36 

412 

23.5 

1.18 

31.08 

503 

23.8 

1.72 

35,11 

1 

596 

26.2 

2.00 

36.73 

635 

27.1 

3.00 

40.87 

738 

29.2 

4.40 

44.33 

'  916 

29.7 

16 


Segmenting  the  splines  was  investigated  to  determine  if  the 
computer  memory  requirements  could  be  reduced  and  to  determine  flutter 
speed  changes.  The  segmented  splines  were  of  two  types:  one  consisting 
of  inboard  and  outboard  sections  (Fig.  16)  and  the  other  consisting  of 
inboard  forward,  inboard  aft,  outboard  forward,  and  outboard  aft 
sections  (Fig.  17).  All  aerodynamic  elements  and  all  structural  nodes 
were  used  as  input  to  each  spline.  It  was  discovered  that  the  flutter 
frequency  and  velocity  did  not  change  with  either  type  of  spline.  The 
amount  of  computer  memory  required  was  reduced  from  300K  for  a  full 
spline  to  200K  for  a  segmented  (inboard  and  outboard)  spline,  therefore, 
the  user  with  limited  computer  resources  could  use  this  procedure  if 
the  splines  are  segmented.  There  was  no  additional  memory  reduction  by 
further  segmenting  the  spline.  Also,  a  full  spline  with  the  number  of 
structural  grid  points  reduced  to  15  (Fig.  18)  was  used  with  result 
showing  a  reduction  in  flutter  speed  and  no  significant  cost  reduction. 

Since  Thomson  (Ref  18)  used  an  aerodynamic  model  consisting  of  20 
aerodynamic  boxes  compared  to  the  64  used  here,  it  was  necessary  to 
determine  the  effect  on  caused  by  such  a  model  change.  With  this 
model  change  and  using  all  structural  nodes,  the  resulting  flutter  con¬ 
dition  was  a  of  846  Knots  and  of  31.6  cps.  This  was  a  significant 
reduction  from  the  916  knot  flutter  speed  previously  discovered.  The 
structural  nodes  input  to  the  spline  were  reduced  to  15  points  and 
using  the  20  box  aero  model,  the  flutter  condition  was  a  of  855  Knots 


and  0)^  was  30.4  cps.  Again,  the  of  855  knots  was  significantly 
below  the  916  knot  flutter  speed.  Based  on  these  results,  the  64  box 
aero  model  is  the  most  appropriate  model  to  use  for  flutter  analysis. 

To  determine  the  effect  on  flutter  speeds,  the  altitudes  were 
varied  and  compared  to  the  sea  level  conditions.  Flutter  will  occur 
at  one  dynamic  pressure  q^.  Since  is  ip„V^  .  the  flutter  dynamic 
pressure  is  considered  a  constant  value.  Therefore,  as  the  altitude 
is  increased,  (density  of  the  free  stream)  will  decrease  and  the 
flutter  velocity  must  increase.  Altitudes  of  20,000  and  40,000  feet 
were  used  in  the  flutter  analysis  procedure  resulting  in  of  1127 
Knots  and  of  30.6  cps  at  20,000  feet.  Flutter  conditions  at  40,000 
feet  were  not  found  using  a  reduced  frequency  range  of  .1  to  1.0. 

These  results  verify  that  does  increase  with  altitude. 

The  range  of  reduced  frequency  (k),  used  in  interpolation  of  the 
aerodynamic  forces,  was  varied  to  determine  if  the  flutter  conditions 
were  sensitive  to  that  parameter.  The  initial  values  of  k  were  select¬ 
ed  such  that  no  interpolation  was  required.  The  reduced  frequencies 
selected  were  decreased  to  two  values  to  span  the  region  of  k  at 
flutter.  This  resulted  in  a  decreased  flutter  speed  with  the  flutter 
frequency  remaining  fairly  constant.  An  increase  in  flutter  speed 
resulted  from  using  three  values  of  k  to  span  a  region  from  .2  to  .9k. 
The  variation  in  flutter  speed  versus  the  range  and  number  of  k  values 
used  are  listed  in  Table  VI. 


Table  VI 


Flutter  Speed  Changes  Due  to  Variation  of  k 


Number  of 
k  Values 


Range  of 
k 


.15  -  .25 
.2  -  .9 
.2  -  .6 


'f  (knots 


“f  (cps 


The  flutter  conditions  obtained  using  5  values  of  k  are  considered  to 
be  the  most  accurate  due  to  the  fact  that  no  interpolation  was  required 
when  computing  the  aerodynamic  forces  at  the  values  of  k  for  the  V-g 
plots. 

All  of  the  factors  affecting  V^  have  been  investigated  resulting 
in  an  aerodynamic  model,  structural  model,  and  splines  with  the  follow¬ 
ing  characteristics:  a  spring  constant  of  4.40  x  10^  in-lbs/rad,  a 
segmented  spline  consisting  of  an  inboard  and  outboard  section  with 
input  points  local  to  each  section,  a  sixty-four  box  aerodynamic  model, 
and  sea  level  conditions.  Using  these  models,  simulated  repairs  and 
their  effects  were  investigated. 
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IV.  FLUTTER  ANALYSIS  OF  A  REPAIRED  MODEL  USING 
THE  DOUBLET  LATTICE  METHOD 

To  properly  model  the  simulated  repairs,  the  T.O.  IT-38A-3  was 
used  to  find  the  current  repair  limitations  and  repair  procedures.  The 
repairs  chosen  for  simulation  were  considered  to  be  the  ones  that  would 
change  the  flutter  characteristics  of  the  stab  the  most.  The  locations 
and  types  of  damage  simulated  were  a  hole  near  the  trailing  edge  (Fig. 
19),  a  hole  at  the  tip  (Fig.  20),  and  delamination  of  the  leading  edge 
where  some  of  the  honeycomb  core  must  be  removed  (Fig.  21).  The  size 
limitations  for  these  areas  (Ref.  19)  are  listed  in  Figures  22  to  25. 
The  size  simulated  repairs  were  modeled  to  be  close  to  the  current 
repair  limitations.  The  repairs  were  modeled  by  changing  the  non- 
structural  mass  of  the  elements  in  the  depicted  repair  areas.  The 
honeycomb  core  has  a  density  of  3.1  lb  per  cubic  foot  and  the  repair 
material  (Nylon  phenolic  core)  has  a  density  of  4.5  lb  per  cubic  foot. 
The  difference  in  the  two  densities  was  added  to  the  repaired  elements 
on  the  structural  model.  Due  to  the  repair  procedures,  no  additional 
mass  was  added  for  the  skin  patch.  Using  the  changed  structural  model, 
the  flutter  analysis  procedure  was  repeated. 

The  in  vacuum  frequencies  of  the  repaired  models  changed  very 
little  from  those  of  the  clean  model.  Since  the  frequencies  showed 
minor  changes  as  a  result  of  the  repairs,  there  was  no  significant 


change  expected  in  the  flutter  conditions.  Table  VII  shows  a  comparison 
of  the  clean  model  to  the  repaired  models. 


With  the  modeled  delamination,  a  slight  increase  in  flutter  speed  was 
obtained.  This  was  due  to  the  mass  addition  forward  of  the  elastic  axis, 
which  has  a  stabilizing  effect.  These  results  show  that  the  flutter 
speed  and  frequency  are  not  greatly  affected  by  repairs  within  the 
current  repair  limitations. 

A  repair,  which  exceeded  the  current  repair  limitations,  was 
modeled  to  investigate  its  changes  to  the  flutter  speed  and  frequency 
(Fig.  26).  This  repair  was  simulated  by  replacing  the  honeycomb  core 


with  the  repair  core  from  the  root  rib  to  the  tip  rib  in  a  section  near 
the  trailing  edge.  The  results  showed  a  flutter  speed  of  911  knots  and 
a  flutter  frequency  of  29.7  cps.  Even  though  this  repair  exceeded 
current  limitations,  there  was  no  significant  change  in  the  flutter 
conditions  from  that  of  a  clean  model. 


V.  FLUTTER  ANALYSIS  USING  STRIP  THEORY 


Utilizing  the  previously  defined  clean  structural  model,  the 
aerodynamic  model  was  changed  to  a  strip  theory  model.  Strip  theory 
is  an  aerodynamic  theory  that  uses  strips  parallel  to  the  free  stream 
direction  to  determine  the  aerodynamic  data  at  subsonic  and  supersonic 
speeds.  The  aerodynamic  model  used  to  represent  the  stab  consists  of 
fifteen  strips  (Fig.  27).  The  spline  must  consist  of  a  representation 
of  all  strips,  therefore,  the  spline  model  and  aerodynamic  model  are 
identical.  All  structural  grid  points  (Fig.  11)  were  used  as  input 
points  to  the  spline.  With  lift  curve  slope  information  along  the  span 
provided  by  NAI  (Ref.  11),  the  flutter  procedure  was  repeated. 

The  flutter  analysis  was  accomplished  at  subsonic  and  supersonic 
speeds.  There  were  no  differences  in  the  flutter  conditions  in  either 
speed  regime,  therefore,  the  results  stated  are  those  computed  at  .5 
Mach.  The  in  vacuum  first  torsion  frequency  of  44.3  cps  was  changed  to 
29.04  cps  at  flutter  with  a  flutter  speed  of  946  knots.  The  structural 
model  was  changed  to  simulate  the  previously  defined  repairs  and  the 
flutter  procedure  was  repeated  and  the  results  are  listed  in  Table  VIII. 


Table  VIII 


Comparison  of  Flutter  Conditions  of  a  Clean  Model 
Versus  Repaired  Models  Using  Strip  Theory 


Repair 

In  Vacuum 

1st  Torsion 

Freq  (cps) 

t 

(knots) 

(cps) 

Clean 

44.3 

946 

29.04 

Hole  T.E. 

44.3 

946 

29.04 

Hole  Tip 

44.3 

946 

29.03 

Delamination 

44.3 

947 

28.98 

Strip 

44.2 

943 

28.98 

0? 


The  results  in  Table  VIII  showed  the  same  trends  found  using  the 
doublet  lattice  method  including  a  slight  increase  in  flutter  speed  for 
the  delamination  at  the  leading  edge  of  the  stab. 


VI.  FLUTTER  ANALYSIS  USING  MACH  BOX 


Another  aerodynamic  theory  considered  was  the  mach  box  method, 
a  supersonic  theory,  which  required  a  change  in  the  aerodynamic  model. 
The  same  clean  structural  model  was  used  in  this  portion  of  the  inves¬ 
tigation.  The  mach  box  aerodynamic  model  consists  of  an  array  of 
rectangular  boxes  whose  diagonals  parallel  the  mach  line  (Ref  1).  At 
the  center  of  each  box  is  a  pulsating  source  emanating  spherical  dis¬ 
turbances.  In  general,  the  model,  defined  by  the  mach  line  and  the 
planform,  consists  of  the  planform,  diaphragm,  and  wake  regions  (Fig. 
28).  In  NASTRAN,  the  user  must  specify  the  number  of  chordwise  boxes 
and  the  mach  number  such  that  the  number  of  boxes  do  not  exceed  200. 

The  aero  model  for  the  stab  using  the  mach  box  method  is  shown  in 
Fig.  29. 

The  spline  must  consist  of  at  least  three  arbitrarily  chosen 
points  on  the  planform.  Six  locations  were  chosen  on  the  stab  to 
represent  the  spline  model.  Inputs  from  all  structural  nodes  (Fig.  11) 
were  used  with  the  six  spline  points  for  use  in  the  interpolation  pro¬ 
cedure. 

Using  sea  level  conditions  and  various  mach  numbers,  the  flutter 
procedure  was  repeated  with  the  results  listed  in  Table  IX. 


Table  IX 


Flutter 

Conditions  Using  Mach 

Box 

Mach  # 

In  Vacuum 

1st  Torsion 
Freq  (cps) 

'^f  (knots) 

‘^U.cp; 

1.10 

44.3 

X 

X 

1.15 

44.3 

X 

X 

1.20 

44.3 

X 

X 

3.0 

44.3 

790 

46.7 

5.0 

44.3 

780 

46.7 

The  results  show  the  flutter  conditions  at  the  higher  mach  numbers  but 
no  flutter  conditions  at  mach  number  1.2  and  below  even  though  the  pre¬ 
dicted  flutter  speed  is  within  that  region.  The  flutter  frequency  of 
46.7  cps  was  not  close  to  the  expected  value  of  29  cps.  Simulated  re¬ 
pairs  were  made  to  the  structural  model  and  the  flutter  procedure  re¬ 
peated.  The  results  show  the  same  trend  discovered  using  the  doublet 
lattice  and  strip  theories,  but,  the  frequency  is  still  far  from  the 
expected  value.  This  data  is  based  on  an  assumed  Mach  number  of  3.0. 


Table  X 


Comparison  of  Flutter  Conditions  of  a  Clean  Model 


Versus  Repaired  Models  Using  Mach  Box 


Repair 
Clean 
Hole  T.E. 
Hole  Tip 
Delamination 
Strip 


In  Vacuum 
1st  Torsion 
Freq  (cps) 


f  (knots! 


46.91 

46.91 

46.91 

46.90 


46.61 


Unlike  the  strip  and  doublet  lattice  theories,  the  results  due  to 
delamination  did  not  show  an  increase  in  the  flutter  speed.  Even 
though  mass  was  added  forward  of  the  elastic  axis,  the  center  of 
pressure  is  at  50  percent  chord  versus  25  percent  chord,  therefore, 
the  result  is  a  decrease  in  flutter  speed.  It  is  the  author's  con¬ 
clusion  that  using  the  mach  box  theory,  there  are  no  flutter  conditions 
in  the  region  of  interest. 


VII.  RESULTS 


The  flutter  speeds  predicted  in  this  report  vary  depending  on  the 
aerodynamic  theory  used.  The  results  from  the  doublet  lattice  method 
are  considered  to  be  the  most  comparable  to  but  still  below  actual. 

This  is  based  upon  a  statement  in  Ref  12  that  the  computed  830  knot 
flutter  speed  is  ten  to  twenty  percent  below  the  expected  value.  Also, 
the  830  knot  flutter  speed  is  based  upon  three  percent  damping  while 
the  computed  flutter  speeds  in  this  report  are  based  upon  no  damping. 
With  the  ten  to  twenty  percent  margin  quoted,  the  flutter  speed  of  the 
stab  should  be  between  913  and  996  knots.  The  916  knot  V^,  predicted 
using  doublet  lattice  method,  is  in  the  lower  portion  of  this  region 
and  is  considered  a  good  estimate  for  the  actual  flutter  speed.  A 
flutter  frequency  of  29  cps  was  reported  in  reference  12,  therefore, 
the  of  29.7  cps  compares  favorably  to  that  value  reported  by  NAI. 

With  the  simulated  repairs  used  in  the  flutter  analysis  procedure,  the 
flutter  conditions  were  virtually  unchanged  from  those  of  a  clean  model. 

The  results  of  this  report  were  obtained  by  using  procedures  that 
deviated  from  those  generally  accepted  for  flutter  analysis.  A  coarse 
aerodynamic  model  is  usually  chosen  but  it  has  been  shown  in  this  in¬ 
vestigation  that  the  flutter  speed  is  sensitive  to  the  size  of  the  aero 
mesh  and  may  change  by  as  much  as  50  knots.  Using  all  of  the 
structural  grid  points  as  input  to  the  spline  versu*-  a  small  number  of 
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grid  points,  the  change  if  was  not  significant  unless  combined  with 
the  changes  in  the  aero  model.  The  best  results  were  obtained  using  a 
fine  aerodynamic  model  and  all  structural  grid  points.  Segmenting  the 
splines  produced  no  change  in  flutter  conditions,  but  resulted  in  a 
significant  decrease  in  the  amount  of  computer  core  memory  required. 

This  reduced  computer  requirement  makes  the  procedures  in  this  report 
available  to  a  larger  number  of  users.  Also,  the  reduced  frequency  is 
usually  chosen  such  that  interpolation  is  required  over  much  of  the 
reduced  frequency  range.  The  flutter  conditions  show  some  inconsis¬ 
tencies  depending  on  the  range  of  k  used,  therefore,  the  k  values  should 
be  chosen  such  that  a  minimal  amount  of  interpolation  for  the  aero 
forces  is  required. 


VIII.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  procedure  used  in  this  investigation  for  flutter  analysis  of 
clean  and  repaired  models  of  the  stab,  ensure  all  factors  affecting  the 
flutter  conditions  are  considered.  According  to  Ref  7,  an  aircraft 
must  be  flutter  free  for  up  to  115  percent  of  the  aircraft  limiting 
airspeed.  For  the  T-38,  the  limiting  airspeed  is  710  knots.  Consider¬ 
ing  the  safety  margin  required  due  to  flutter,  the  aircraft  should  be 
flutter  free  up  to  an  airspeed  of  816  knots.  From  the  results  of  this 
investigation,  the  changes  in  due  to  the  repairs  are  shown  to  be 
minor;  therefore,  it  is  recommended  that  the  procedures  in  this  report 
be  used  to  determine  the  possibility  of  establishing  new  repair  limit¬ 
ations  for  the  stab.  An  additional  safety  factor  can  be  added  to  the 
916  knot  speed  to  ensure  that  speed  is  not  closely  approached. 

It  is  recommended  that  the  flutter  procedure  be  initiated  using 
a  coarse  k  set  in  order  to  determine  the  reduced  frequency  region  at 
which  flutter  will  occur.  Since  it  has  been  shown  that  the  flutter 
speed  varies  depending  on  the  reduced  frequency  used,  the  coarse  k  set 
should  be  modified  to  a  region  in  the  vicinity  of  the  expected  reduced 
frequency  at  flutter.  This  modified  k  set  should  be  chosen  such  that 
minimal  interpolation  for  computation  of  the  aerodynamic  forces  is 
required. 


For  establishing  new  repair  limitations,  it  is  recommended  that 
the  doublet  lattice  method  be  used.  Within  the  current  repair  limita¬ 
tions,  it  appears  as  though  strip  theory  gives  the  best  results  but  as 
more  mass  is  added  near  the  trailing  edge,  there  should  be  chordwise 
deformation  and  as  a  result,  the  flutter  speed  is  affected.  Chordwise 
deformation  is  not  detected  by  strip  theory,  therefore,  there  is  a 
larger  decrease  in  flutter  speed  using  strip  theory  when  compared  with 
the  flutter  speed  decrease  using  the  doublet  lattice  method. 
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Figure  2  NASTRAN  T-38  Series  3  Stabilizer  Model 


Figure  3  NASTRAN's  1st  Bending  Mode,  17.36  cps 


Figure  5  NASTRAN's  1st  Torsion  Mode,  44.42  cps 


Figure  6  NASTRAN's  2nd  Bending  Mode,  77.58  cps 


3rd  Bending  Mode,  124.68  cps 


Figu»-e  9  NASTRAN's  3rd  Torsion  Mode,  170.38  cps 
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Figure  13  NAI' 


Figure  14  V-6  Diagram  of  NASTRAN  Computer  Modes  4,5,6 


■e  23  Repair  Limit,  Hole  Tip  (Ref  19) 
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e  29  Mach  Box  Aerodynamic  Model  of  the  Stabilizer 
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